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ABSTRACT 

In  a  liquid  a  satisfactory  method  of  defining  the  diffusivity  so  that  it  will  be  independent  of  relative  composition, 
has  yet  to  be  devised.  The  recent  development  of  effective  methods  for  the  calculation  of  mass  -transfer  rates  from  particle 
surfaces,  now  brings  forward  the  need  of  considering  pore  diameters,  permeabilities  and  porosities.  Knowledge  of  liquid 
diffusion  rates  in  porous  solids  has  applications  in  solid-liquid  extraction  and  in  catalysis.  Piret  E.  L.,  Ebel  R.A.  Kiang  C.T. 
Armstrong  W.  P.  aimed  to  investigate  about  liquid  diffusion  rates  which  would  be  of  value  in  the  study  of  actual  systems 
involving  porous  carrier  solids  and  to  obtain  idealized  extraction  systems.  In  their  work,  they  presented  theoretical 
expressions  and  experimental  data  for  extraction  from  single  capillaries,  from  beds  of  uniform  glass  beads  and  from  inert 
porous  spheres  carrying  single-phase  solutions  of  soluble  substances.  They  introduced  the  concept  of  pore  shape  factor. 
Pore  shape  factor  is  equal  to  the  ratio  of  the  actual  to  the  theoretical  extraction  time  or  for  an  actual  sphere  of  measured 
radius  R,  an  equivalent,  idealized  sphere  of  radius  KR,  (K  >  1.0)  can  be  postulated,  whose  structure  offers  no  resistance  to 
diffusion.  The  square  of  this  constant,  K2,  will  be  called  the  pore  shape  factor.  This  factor  is  used  in  interpreting  and 
correlating  the  extraction  data.  Here  in  this  paper  the  author  intends  to  interpret  the  results,  obtained  in  their  investigations 
and  to  derive  a  simple  formula  to  determine  this  pore  shape  factor.  For  this  purpose,  Euler's  Summation  Formula  and 
Abel's  identity  are  applied  to  the  analytic  solutions  obtained  for  their  experiments.  This  approximation  is  obtained  only  for 

small  values  of  -j,  where  D  is  diffusivity,  t  is  time  and  R  is  measured  radius  of  sphere.  The  formula  is  verified  for 

experimental  data.  It  shows  nice  agreement  with  the  experimental  data. 

KEYWORDS:  Extraction,  Pore  Shape  Factor,  Euler's  Summation  Formula,  Abel's  Identity  and  Big  oh  Notation 
INTRODUCTION 

Related  to  investigations  of  extraction  rates,  a  series  of  experiments  was  carried  out  on  the  batch  extraction  of 
single  phase  solutes  from  porous  spheres  by  Edgar  L.  Piret,  R.  A.  Ebel,  C.T.  Kiang  and  W.P.  Armstrong  (University  of 
Minnesota,  Minneapolis,  Minnesota).  [1]  Aqueous  solutions  of  potassium  chloride,  sodium  chloride,  potassium  chromate, 
copper  sulfate,  and  acetic  acid  were  extracted  from  three  types  of  carrier  solids  of  increasing  complexity:  single  capillaries, 
tubes  containing  beds  of  small  glass-beads  and  porous  alumina  spheres.  The  transfer  mechanism  was  verified  as  being 
diffusional.  Equations  were  derived  for  the  extraction  of  a  single  solute  from  single  capillaries  and  from  spheres  in  the 
batch-wise  operation.  Validity  of  the  equation  of  batch-wise  operation  was  verified  by  extractions  from  porous  alumina 
spheres  under  several  conditions  of  concentration,  solvent  volume  and  temperature. 

The  influence  of  an  inert  porous  solid  upon  diffusion  rate  has  not  been  extensively  investigated.  Cady  and 
Williams  [2]  studied  diffusion  of  urea  through  wood  fibers.  They  stated  that  the  diffusion  rate  was  decreased  principally  by 
a  mechanical  blocking  effect  of  the  cell  structure;  the  effect  of  cell  wall  friction  on  molecular  motion  was  of  much  less 
importance. 
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Boucher,  Brier  and  Osburn  [3]  found  that  the  extraction  of  soybean  oil  from  porous  chaplets  was  a  purely 
diffusional  operation  whose  rate  was  decreased  by  the  nature  of  the  solid.  They  suggested  the  factor  affecting  rate  to  be  the 
mechanical  blocking  by  the  solid,  the  drag  on  molecule  due  to  the  proximity  to  cell  walls  and  the  tortuous  path  through  the 
solid. 

Piret  E.  L.,  Ebel  R.A.  Kiang  C.T.  Armstrong  W.  P.  studied  several  facts,  which  affect  liquid  diffusion  rates  and 
which  would  be  of  value  in  the  study  of  actual  systems  involving  porous  carrier  solids.  In  their  experiments,  they  studied 
extraction  from  single  capillaries,  from  beds  of  uniform  glass  beads  and  from  inert  porous  spheres  carrying  single-phase 
solutions  of  soluble  substances.  The  new  thing  is  that,  they  introduced  the  concept  of  pore  shape  factor.  Pore  shape  factor  is 
equal  to  the  ratio  of  the  actual  to  the  theoretical  extraction  time  or  for  an  actual  sphere  of  measured  radius  R,  an  equivalent, 
idealized  sphere  of  radius  KR,  (K  >  1.0)  can  be  postulated,  whose  structure  offers  no  resistance  to  diffusion.  The  square 
of  this  constant,  K2,  will  be  called  the  pore  shape  factor.  This  factor  is  used  in  interpreting  and  correlating  the  extraction 
data. 

RESEARCH  OBJECTIVES 

It  was  the  purpose  of  the  work  of  the  team  of  Piret  E.  L.,  Ebel  R.A.  Kiang  C.T.  and  Armstrong  W.  P.  to  obtain  on 
idealized  extraction  systems  further  basic  information  on  liquid  diffusion  rates.  They  conducted  experiments  using:  (1) 
Single  Capillary  Investigations  and  (2)  Porous  Sphere  Investigations. 

In  this  paper  the  data  obtained  during  their  investigations  are  interpreted  mathematically  and  a  simple  formula  to 
predict  pore  shape  factor  is  derived. 

SINGLE  CAPILLARY  INVESTIGATION 
Procedure  and  Data 

A  straight  capillary  tube  of  constant  bore  and  known  internal  length  and  volume,  but  with  one  end  sealed  off,  was 
filled  with  a  binary  solution  of  known  composition.  The  capillary  was  held  in  a  vertical  position  with  the  open  end  up,  and 
a  very  slow  stream  of  pure  solvent  was  allowed  to  flow  at  right  angles  to  the  longitudinal  axis  of  the  capillary.  The 
temperature  was  maintained  at  25.0  ±  0.1  °C.  After  an  elapsed  time  t,  the  capillary  was  removed  and  the  solution  inside  the 
capillary  was  washed  into  a  flask  or  crucible  and  analyzed  volumetrically  or  gravimetrically.  Amount  and  composition  of 
the  solution  now  inside  the  capillary  were  determined.  Generally  several  such  experiments  were  performed  and  the  data 
presented  as  E  (  per  cent  of  extractable  material  remaining  unextracted  )  vs.  time  of  extraction  t. 

They  used  the  capillary-tube  method  to  determine  the  integral  diffusivity  for  a  number  of  systems.  Their  data  for 
KCl  —  H20  solutions  (C0  =  4  N)  are  given  below: 


Table  1 


No. 

ti      sec / 
l\2  '  /cm2 

( 100  -  £)% 

1 

51.87  x  103 

95.1 

2 

36.50 

85.8 

3 

28.97 

80.4 

4 

17.03 

63.2 

5 

8.805 

45.9 

6 

3.689 

28.8 

7 

2.838 

24.9 

8 

1.242 

15.9 

9 

1.160 

15.8 

10 

0.741 

12.2 
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Interpretations 

Assuming  that  the  solvent  concentration  is  constant,  and  CA  is  the  solute  concentration  and  D  is  diffusivity  of  a 
binary  liquid  system,  the  governing  partial  differential  equation  will  be, 

d_cA=^rDdcA 

dt        dx\     dx  J  w 

with  the  boundary  conditions 
CA  —  C0  at  t  —  0,0  =  x  —  L 
CA—  0  at  x  —  L,t  >  0 

Dir=hftc*dx;  atx  =  L>t>0 

Assuming,  diffusivity  D,  to  be  constant,  the  solution  [1],[4],[5]  to  equation  (1),  will  be 

800  Y30  1  (-{2n-l)Vpt\ 

E  -  -rLn=i^-^exp  \  —2  )  (2) 

For  sufficiently  large  values  of  t,  equation(2)  may  be  approximated  as 

(Considering  only  first  two  terms,  as  contribution  from  the  terms  after  that  will  be  negligible.) 
Taking,  y  —  exp  (  "Lft),  equation  (3)  results  in,  E  —  ^  [y  +  ^y9 j 

Applying  this  equation  for  the  data  for  which  E  is  between  25%  to  75%,  That  is  for: 
(100  -E)  =  63.2, 45.9,  28.8,  (for  fourth,  fifth  and  sixth  observation  in  table), 
y  =  0.4539, 0.6648, 0.8488,  respectively. 

/-n2Dt\  4xlny 

Moreover,  y  =  exp  [—)  =>  D  =  . 

So,  diffusivity  D  =  1.87977  X  10-5, 1.87922  X  10"5, 1.801  X  10"5  cm2 /sec  respectively. 

Based  on  the  capillary  experiments,  they  determined  integral  diffusivity  of  potassium  chloride  to  be  1.87  X 
10"5  cm2 1  sec.  (Literature  value  of  diffusivity  of  KCl  is  1.84  X  10"5  cm2  /sec. ) 

POROUS  SPHERE  INVESTIGATIONS 
Procedure  and  Data 

A  series  of  experiments  was  carried  out  on  the  batch  extraction  of  single  phase  solutes  from  porous  aluminum 
spheres.  This  was  done  by  filling  capillaries  of  the  porous  spheres  by  salt  solutions  and  determining  the  rate  of  extraction 
into  a  measured  volume  of  solvent.  The  purpose  of  these  experiments  was  to  determine  the  effect  of  porous  medium  upon 
extractions.  Single-phase  solute  is  contained  in  porous  sphere  and  the  sphere  is  placed  into  a  finite  quantity  of  agitated 
solvent.  As  the  diffusion  proceeds,  the  concentration  in  the  outside  solution  rises. 

The  following  table  shows  their  experimental  data  for  extraction  of  aqueous  Potassium  Chloride  from  porous 

spheres: 
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Table  2 


Salt  KC1 

Initial  Outside  Concentration  =  0 

Temperature"  c 

25 

Initial  Inside 
Concentration  g/ml 

0.114 

25 

t  (Minutes) 

(100 

-E)  % 

t  ( Minutes) 

(100-E)  % 

t  ( Minutes) 

(100-E)  % 

5 

16.3 

10 

24.1 

20 

33.8 

15 

28.2 

30 

40.4 

60 

49.4 

45 

46.9 

60 

54.0 

150 

71.6 

240 

82.6 

180 

78.0 

240 

82.3 

300 

85.4 

375 

90.8 

360 

89.5 

480 

92.3 

540 

95.6 

600 

95.0 

1140 

98.8 

1380 

99.0 

720 

95.4 

Interpretation 

The  equation  for  symmetrical  diffusion  from  a  porous  sphere  of  porosity  P  is 


That  is 

dc       „  (d2c      2  dc\ 


dPc  _  fd2Pc  2  dPc 
dt        V  dr2     r  dr 


Kdr2      r  dr) 


=  d(° 

dt  \d 

Defining  a  new  variable, 
u  =  cr,  (4)  converts  into 

dt  \dr2J 

The  boundary  conditions  are, 
u  =  0  at  r  =  0 
u  =  c0r  when  t  =  0 
u  =  cr  when  t  >  0 
c  =  cx  =  Co,  when  t  =  oo 

where  c  —  c(  r,  t)  =  Concentration  of  inside  solution  at  any  point  r  and  time  t ,  g./cc 

C  —  Concentration  of  outside  solution  at  any  time,  g./cc 

R  =  Radius  of  the  sphere 

The  solution  [1],[5]  to  this  problem  is 

f-mp2Dt\ 


_  200o2 
£  -       3  2wp=l 


(4) 


(5) 


/  amv2\ 


(6) 


Introduction  to  Concept  of  Pore-Shape  Factor  for  Porous  Spheres  in  the  Study  of  Extraction  Rates 


97 


Here  mp '  s,  are  positive  roots  of  the  equation  tanx  = 


X 


a  —  solution  ratio,  volume  of  outside  solution  per  unit  volume  of  solution  in  porous  sphere. 
K2  =  pore  shape  factor,  to  be  determined  yet. 

In  practice,  porous  spheres  with  all  these  characteristics  were  difficult  to  be  obtained.  Here  an  actual  sphere  of 
radius  R  which  usually  had  tortuous,  constricted  capillaries  is  considered  to  be  equivalent  to  an  idealized  sphere  of  radius 
KR.  (K  >  1.0)  The  increased  length  of  path  and  the  throttling  effect  of  constrictions  provide  resistance  to  diffusion.  In 
fluid  flow  a  valve  is  considered  to  be  equivalent  in  resistance  to  a  certain  length  of  straight  pipe.  Here  an  actual  sphere  of 
radius  R  which  usually  has  tortuous,  constricted  capillaries  is  considered  to  be  equivalent  to  an  idealized  sphere  of  radius 
KR.  (  K  >  1.0  )  The  square  of  this  constant,  K2  will  be  called  the  pore  shape  factor.  This  factor  is  also  equal  to  the  ratio  of 
the  actual  to  theoretical  time,  based  on  an  ideal  sphere,  for  a  given  per  cent  extraction. 

In  their  experiments,  they  made  early  runs  with  an  initial  inside  concentration  of  zero  and  a  3Af  outside 
concentration  for  the  purpose  of  determining  pore  shape  factors  of  alumina  spheres.  Results  were  plotted  on  semi-log  paper 
with  E  per  cent  unextracted  as  the  ordinate  and  time  t  in  minutes  as  the  abscissa.  They  drew  curves,  taking  K2  —  4,  K2  — 
7,  K2  —  9,K2  —  1A  ,  and  substituting  an  integral  diffusivity  of  KCl  to  be  1.87  X  10"5  sq.  cm./sec  ,  which  was 
determined  by  capillary  experiments,  in  equation  (6).  For  large  value  of  t,  experimental  results  were  close  to  curve  for 
K2  =  9,  otherwise  all  points  fell  close  to  a  single  theoretical  curve  for 

K2  =7.  [1]  Thus  it  may  be  seen  that  variations  of  inside  concentration  have  no  effect  on  the  data.  This  is  an 
agreement  with  equation  (6).  They  accepted  value  of  K2  to  be  7. 

Results  may  be  verified  in  other  way  also  as  follows:  For  large  value  of  t,  equation  (6)  can  be  approximated  as, 


considering  only  first  term  of  infinite  series,  because  then  onwards  contribution  from  the  terms  will  be  negligible. 


(7) 


For  a  =  25,  first  positive  root  of  tanx  = 


i+ 


~25~z  is  nix  =  3.179  (It  is  close  to  ) 

 X2 


3 


Substituting  a  =  25,  D  =  1.87  x  10~5  cm2 


I  sec  ,R  =  0.5985  cm,  (as  given  in  [1]) 


and  m±  —  3.179  in  (7),  now  expression  will  be 


(8) 


Using  (8),  the  obtained  results,  are  as  follows. 
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Table  3 


t  (in  Minutes) 

t  (in  Seconds ) 

E  % 

K2 

240 

240x  60 

17.4 

6.379 

300 

300x  60 

14.6 

6.950 

375 

375x  60 

9.2 

6.493 

480 

480x60 

7.7 

7.574 

1140 

1140x60 

1.2 

9.337 

It  can  be  noticed  that  the  values  of  K2  obtained  here,  are  close  to  7,  for  the  first  four  values  in  the  table  and  close 
to  9,  for  the  last  value,  similar  to  the  results,  obtained  using  theoretical  curves. 

DERIVATIVE  OF  A  SIMPLE  FORMULA  TO  DETERMINE  K2 

Precession  and  accuracy  seem  better  for  small  values  of  t.  Instead  of  using  approximation  for  large  value  of  t, 
approximation  for  small  values  of  t,  would  be  better.  Here  in  this  paper  the  author  intends  to  derive  a  simpler  formula  to 
determine  value  of  K2.  This  formula  will  be  valid  for  only  small  values  of  ^Vrz-  ^  w^  derived  taking  help  of  Euler's 
Summation  Formula  and  Abel's  Identity.  Statements  of  these  theorems  [6]  are  as  follows: 

Euler's  Summation  Formula 

If  /  has  a  continuous  derivative  /  on  the  interval  [  y,  x  ]  ,  where  0  <  y  <  x  ,  then 


Zy<„s*/(n)  =  S'fWt  +  /*(  t  -  [t])fXt)dt  +  /(*)([*]  -  x)  -  f(y)(]y]  -  y) 
Abel's  Identity 

For  any  arithmetical  function  a(n),  let  A(x)  =  Y*n<x  a(n),  where  A(x)  =  0,  if  x  <  1. 
Assume  /  has  a  continuous  derivative  on  the  interval  [  y,  x],  where  0  <  y  <  x.  Then 

Tjy<n  <x  1 

Derivation 


(9) 


.a(n)f(n)  =  A{x)f{x)-A{y)f{y)  -  f* A(t)f\t)dt. 


(10) 


With  the  aim  to  apply  these  wonderful  theorems  to  the  solution  (6),  here  it  is  required  to  determine  an  arithmetical 
function,  which  can  be  utilized  as  a(n)  and  a  function  /  , which  has  continuous  derivative.  For  a  =  25,  values  of  first  10 


roots  of  tan  x  — 


are: 


m1  =  3.179,m2  =  6.302,  m3  =  9.438, m4  =  12.576,  m5  =  15.715, 

m6  =  18.856,m7  =  21.997,  m8  =  25.137,m9  =  28.275,m10  =  31.420. 

It  can  be  noticed  that  as  n  becomes  large,  mn  -»  nn.  Actually  the  difference  between  m2  and  2n  is  also  less  than 
0.018815  and  gradually  the  difference  goes  on  decreasing,  remaining  positive,  as  n  increases. 


So, 


can  be  approximated  by 


Ha(l+^)J 


That  is,  by 


(l+a)+^-n27T2 
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1  19  19  7Z2  3 

y  t  ;  2\T  —  £  ~~2         =  "T^Hi0^  =  ~^r^  x  —  =  — 7,  confirms  its  convergence. 

1 1  a(l  I amn  J  — n2n2         a  n         n         a  n        6  2a 

So,  here      — can  serve  as  a(n).  As  exponential  functions  have  continuous  derivatives, 


exp  (  nK2R2t')  can  serve  nere  as  /(n)- 


Let  y  be  a  number  between  0  and  1.  So,  A(y)  =  £n<ya(n)  =  0,  which  fulfills  requirement  demanded  in  the 
statement  of  the  identity. 

Let  x  be  a  very  large  number.  So 

AM  =  Enfi*a(n)  =  X=1a(n)  -  ^>xa(n) 

_  2a2        a27r2  ^"^n2 

Now  according  to  Euler's  summation  formula 

•5-,  1         rM  dt  rMt-\t\    .      .    [M]-M       \x]-x       ....  ,  ,  ,.  .    „  . 

Lx<n<M  ^  =  ix  ^~2)x  —dt  +  —2  —>  where  M  1S  verY  larSe  number,  tending  to  infinity. 


=  _  (  ^  _  ~)  +  0(x"2)  ~  ^  as  M  is  tending  to  infinity. 
Here-2/M^dt  +  M^l_M^  =  o(x"2), 

Jx     t3  M2  x2  v  J 


means  that  —  -JC  £L_  is  bounded  for  all  x  >  0. 


The  symbol  0(x  2)  represents  an  unspecified  function  of  x  which  grows  no  faster  than  some  constant  times  x  2 . 
The  symbol  N0  'is  known  as  big  oh  notation  [6]. 

Thus  40)  =       -  -^(-)  (11) 

v  J       2a2        a2n2  \xj  v  ' 

Now  using  (11) 

„  9  /  n2jr2Dt\ 

2,y<n<x  a2n2jr2  e*P  ^-  ~^TJ 

(x27T2  Dt\  /     y2n2  Dt\ 

— k2r2  )  ~*  ^'^S  ^  *s  ver^  sma^  number  exP  I — ^2^2  )  ->  1. 

Here  4(p)  =  ^  -  ^  g)  (using  (11))  and  /  '(p)  =  exp  (-  ^)  X  ^g^. 

On  substituting  these  results,  now  (12)  yields, 

Z9  /    n2n2Dt\  fx  (  3  9  /    p2n2Dt\  -2pn2Dt 

y<n<x  '  y    \  / 

-  _  f*f—  exv  {-  p2n2m)  x  -2pn2°^dv  +  (x  —  (  exv  (-  ^E£\  v  "^^^  dv 

~        )y  l2a2    eXP  {         K2R2   )  X      K2R2     )«P    +    )y  a2n2  \  eXP  {         R2R2    )  X     R2R2    )  W 
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As  x  is  very  large  and  y  is  very  small,  there  is  no  harm,  if  x  is  replaced  by  oo  and  y  is  replaced  by  0+.  Hence  (13) 


yields 


Z_i  a2n2n 


n=l 


9  /  n2n2Dt 


\  3  /    p2n2Dt\  9     f°°       /  p27 


p2n2Dt\  -2n2Dt 


i-f  x^(2Jo;exp(-s^s) 


3        Vot w  9 

2  a2  ATfl 


x  V^as  (2  JQ+ exp(-s2)ds)  =  r(1/2)  =  Vtt 


3  9  DC 

2a2  a2K^JnR2 


(14) 


Using  (14)  in  (6),  result  will  be  £  = 


200a^  /  3 


9         DC  \         .  _  .       600  Dt 
3      \  2a2       a2K^nR2  K  -Jjrfl2 


That  is,  £  =  100  - 


600  /  Dt 
if  AjJTfl2 


(15) 


As  the  minimum  possible  value  of  E  is  zero,  the  value  of  —  /-^  can't  exceed  100, 


so  this  formula  is  valid  for  only  small  values  of 


IDC 

R2 


,  600       Dt       ,    .  I  Dt       .   K  Dt     .  nK 

and —  — r  <  100  =>  — r  <  -  =>  —  <   

K   AjJTfl2  AjJTfl2  6  R2  36 


/o7       7         (600)2XDt/  2 

Thus,  for  small  values  of  \—,KL—   — — 

Ajfl2  (100-£)2X7T 


Using  (17)  for  small  values  of  t  (  t  <  60  )  frm  Table-2,  results  obtained  are  as  follows. 

Table  4 


Initial 

t(in 

Dt 

(100 

—  £)% 

K2 

(1)  0.224 

5  x  60 

0.0156615 

16.3 

6.75479 

15x  60 

0.0469846 

28.2 

6.77033 

45x  60 

0.1409539 

46.9 

7.34318 

(2)  0.114 

10X  60 

0.0313231 

24.1 

6.17992 

30x  60 

0.0939693 

40.4 

6.59744 

60x  60 

0.1879385 

54.0 

7.38551 

(3)  0.307 

20x  60 

0.0626462 

33.8 

6.28368 

60x  60 

0.1879385 

49.4 

8.82499 

(16) 
(17) 


Introduction  to  Concept  of  Pore-Shape  Factor  for  Porous  Spheres  in  the  Study  of  Extraction  Rates  101 

CONCLUSIONS 

•  It  can  be  noticed  that  except  for  last  observation  in  Table-4,  all  the  experimental  data  show  good  agreement. 
There  might  be  any  error  in  recording  the  last  observation. 

•  Taking  the  mean  of  first  seven  observations,  K2  =  6.75926  with  standard  deviation  s  —  0.46926. 

•  It  can  be  also  noticed  that  initial  inside  concentration  has  no  effect  on  extraction  rate. 

•  It  can  be  noticed  that  if  ^|  <  0.09,  that  is  if  t  <  30  minutes,  then  there  is  better  agreement  in  results. 

n  Dt 

Thus  Formula  (17)  can  be  used  to  predict  pore  shape  factor  K  for  small  values  of  — . 
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